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We derive analytic formulae for the energy average (including the energy average of the fluctuation 
contribution) and variance of the intraband decay intensity of a superdeformed band. Our results 
may be expressed in terms of three dimensionless variables: F^/Fs, Fjv/d and Fj^/(Fs +1^). Here 
is the spreading width for the mixing of a superdeformed (SD) state |0) with the normally 
deformed (ND) states \Q) whose spin is the same as |0)'s. The \Q) have mean lever spacing d and 
mean electromagnetic decay width Fjv whilst |0) has electromagnetic decay width Fs- 

The average decay intensity may be expressed solely in terms of the variables F^/Fs and F^/d 
or, analogously to statistical nuclear reaction theory, in terms of the transmission coefficients To (15) 
and Tv describing transmission from the \Q) to the SD band via |0) and to lower ND states. 

The variance of the decay intensity, in analogy with Ericson's theory of cross section fluctuations, 
depends on an additional variable, the correlation length Fn/(Fs + r ) = ^Tn/(Fs + F^). This 
suggests that analysis of an experimentally determined variance could yield the mean level spacing 
d as does analysis of the cross section autocorrelation function in compound nucleus reactions. 

We compare our results with those of Gu and Weidenmiiller. 



I. INTRODUCTION 

The basic feature to be explained in the decay-out of superdeformed (SD) rotational bands jl], ^, |[ || is that the 
intensity of the collective 7-rays emitted during the cascade down an SD band remains constant until a certain spin 
is reached where-after it drops to zero within a few transitions. The sharp drop in intensity is believed to arise from 
mixing of the SD states with normally deformed (ND) states of identical spin §. The model of Refs. jf, @, |, § 
attributes the suddenness of the decay-out to the spin dependence of the barrier separating the SD and ND minima 
of the deformation potential. Refs. |To[ O, [l2| discuss the effect of the chaoticity of the ND states on the decay-out. 

In the present paper we derive analytic formulae for the energy average (including the energy average of the 
fluctuation contribution) and variance of the intraband decay intensity of a superdeformed band in terms of variables 
which usefully describe the decay-out. || [|, [l3|, [l4|. We achieve this using the MIT approach to statistical nuclear 
reaction theory ||, |f, fL7j. The MIT approach tackles calculation of the fluctuation cross section and other moments 
of the 5-matrix by directly calculating averages of fluctuating functions of energy. 

In agreement with Gu and Weidenmiiller [jl3| (GW) we find that average of the total intraband decay intensity 
can be written as a function of the dimensionless variables T^/Ts and T^/d where is the spreading width for the 
mixing of an SD state with normally ND states of the same spin, d is the mean level spacing of the latter and T$ 
(Tjv) are the electromagnetic decay widths of the SD (ND) states. Our formula for the variance of the total intraband 
decay intensity, in addition to the two dimensionless variables just mentioned, depends on the dimensionless variable 
IV/fTs + r^)- This additional variable is analogous to the correlation width of the Ericson's theory of cross section 
fluctuations. Its appearance suggests that measurement of the variance of the decay intensity could yield the mean 
level spacing d. 

The paper is organised in the following manner. In Section Q| we express the average decay intensity as an average 
term plus the average of a fluctuation term. In Section |ni] the expression for the average decay intensity obtained 
in Section [ij is evaluated approximately. In Section 1^ we calculate the variance of the intensity within the same 
approximation scheme as Section III. In Section W we interpret our results by analogy with statistical nuclear 
reaction theory, expressing our results in terms of transmission coefficients and compare our results with those of 
GW. In particular we suggest a reason why GW did not observe the variable rAr/(r,g + T*) in their calculation of the 



variance. Finally in Section VI we make some concluding remarks including mention of the limitations of our results. 



II. ENERGY AVERAGES AND FLUCTUATIONS 

The total intraband decay intensity has the form O, O, n8[ 



/oo 
dE\A 00 (E)\ 2 , (1) 
-OO 
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where the intraband decay amplitude is given by 

Aoo(E) =7s<0|G(£)|0) 7 s. (2) 

Here 75 is the electromagnetic decay amplitude of superdeformed state |0) defined such that Ts — Jg- We assume 
that the amplitude feeding |0) is also given by 7s. The total Green's function, G, is given by 

G{E) = (E-H)-\ (3) 

where the full nuclear Hamiltonian is denoted by H and has a non-Hermitian part which accounts for coupling to the 
electromagnetic field [see Eqs. (|), (|l8|) and © below]. The tilde is used to indicate the dual state or adjoint Jl^] of 
|0). 

In what follows we employ the optical background representation introduced Kawai, Kerman and McVoy [fL6|. 
These authors investigated fluctuation cross sections using a representation of the S-matrix in which the background 
S-matrix is chosen to be the energy average of the S'-matrix itself, that is, the S'-matrix corresponding to the optical 
potential. Here we use the same idea to decompose thedecay amplitude, Eq. (^|), into the sum of its Lorentzian energy 
average (Lorentzian energy averaging interval /) , 



A 00 {E) = A 00 {E + iI/2), (4) 

plus a fluctuating part. 

We proceed by introducing Feshbach's projection operators [jl5| 

P = |0)(0| and Q =1— P. (5) 

Let us introduce the notation Q = PGP for the effective Green's function in the P-space and let Hpp = PHP, 
Hpq = PHQ etc. Using the techniques of Ref. |L6| we obtain 

g = g + g a , (6) 

where the average effective Green's function Q = Q(E + H/2) is given by 

G=[E-H\~ 1 (7) 

and the average effective Hamiltonian is given by 

H = Hpp+Wpp, (8) 
where coupling to the eliminated Q-space is accounted for by 

1 _2i 

Wpp = Hpq e-h qq + u/2 Hqp = — Vp Q y Q p ^ 



and the energy dependent coupling potential, V, introduced in Ref. |16 is defined by 

VpQ = H P Q\ 



U/2 



E-Hqq+U/V 



Vqp - V E-Hqq + 1 I/2 Hqp - 
The fluctuating part of the effective Green's function, Q, is given by 

g tt = gv P o Vo P g, (io) 

y E — Hqq — Wqq QPy ' ^ ' 

where coupling back to the P-space is accounted for by 

w QQ = v QP gvpQ. (ii) 

By construction C? fl — in so much as g is unchanged by reaveraging, that is if g = g. 
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The corresponding decomposition of Eq. (g) for the transition amplitude is 

A 00 (E) = A^(E) + A% (E), (12) 
where the energy average of the transition amplitude is 

Aon = isGools (13) 
and the fluctuating part of the transition amplitude is 

4>o=7s0oo7s. ( 14 ) 
Thus the average of Eq. (|lj) for the relative intensity may be written as the incoherent sum 

-fin = Itn + lfL> (15) 

where 

, f°° 2 

C = = (^Ts)- 1 / dE Aoo(E) (16) 

./ — OO 

and 



Jf n = (2 7 rr s )- 1 / dE\AUE)\\ (17) 

<J — OO 

Up to this point no assumptions have been made except that the transition amplitude can be written in the form of 



Eq. (|2j). As will be made clear in Section 111 the meaning of this assumption is that |0) is a doorway for the decay from 
the SD band to the ND states and vice versa. In the manipulations subsequent to Eq. (|J) we have put the average 
intensity in a form consisting of a background term coming from the smooth energy dependence of the doorway plus 
a term resulting from fluctuations on this background. 

The representation we have used is particular suitable for approximation when the ND states are overlapping. In 



Section III we evaluate Eqs. (|15|-|l7|) for the average decay intensity, 7; n , assuming that this is the case. In Section IV 



we calculate the variance which describes the way in which 7j n fluctuates about I m . 

III. AVERAGE DECAY INTENSITY 

Let us assume that Hpp satisfies the eigenvalue equation 

H PP \0) = (E -iTs/2)\0) (18) 

and Hqq 

H QQ \Q) = {E Q -iT N /2)\Q). (19) 

Here, Eq denotes the energy of SD state |0), T$ it's electromagnetic width for decay to the next lowest state in the 
SD band, Eq (Q = 1,...,N) the energy of the N ND states \Q) with the same spin as |0) and Tjv the common 
electromagnetic width of the \Q) for decay to ND states of lower spin. 
Further, let us write the matrix element of Wpp, Eq. (|9[), as 

W 00 = (0\Wpp\0) = A 1 -iT L /2. (20) 

Here, A 1 = ReTFoo is an energy shift which we ignore and 1^ = — 2ImI / l /r oo- Combining these definitions with Eqs. (Q), 
(|J) and (|l^) the average of the transition amplitude can be written 

^ = E-Eo + Ws + TW (21) 

We see that Eq. ( pi] ) exhibits the structure of an isolated doorway resonance. The doorway |0) has an escape width 
Ts for decay to the SD state with next lower spin and a spreading width for decay to the ND states with the same 
spin which are reached by tunnelling through the barrier separating the SD and ND wells. The doorway structure 
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of Eq. ( |21| ) is due to the assumption that the transition amplitude can be written as in Eq. (|2|). The most general 
expression for the transition amplitude has the form A a b — Jab + J2cc' lac{c\G\c')^ c n, where 7 ac describes the coupling 
of channels a and c. In our doorway model the c and c' stand for |0) or \Q), Q = 1, ...,N, ie. ^ c l c )( c l = P + Q- an d 
a and b denote channels the (electromagnetic) coupling to which is taken is accounted for by the non-Hcrmitian part 
of H [Eqs. ( |l8| ) and dTa)], that is they denote the SD state above |0), the SD below |0) and ND states whose spin is 
different from that of |0). The direct coupling of channels a and b, jab, is taken to be zero. 

In order to evaluate Aq it is useful to introduce eigenvectors and eigenvalues of the operator Hqq + Wqq defined 

by 

(H QQ + W QQ )\q) = (E q -iT q /2)\q), q=l,...,N. (22) 
Then from Eq. @ and Eq. © for G R we obtain 

4 fl ~ \- (0\VpQ\q)(q\VQp\0) r „ 

A 0Q - 1S Q^ E -E q + l T q /2 GooJS 

-^E g g t^/^ (23) 

9o q = (0\V PQ \q), (24) 
g q o = (q[V QP \Q). (25) 

We now employ some statistical assumptions which are frequently used in statistical nuclear reaction theory |l5| to 
derive an analytic formula for the decay intensity I m . We shall assume that the Lorentzian and box energy averages 
and the average over the label, q, are all approximately equal, ie., that for a suitable function f q (E) of q and E 



where 



i r E a +AE/2 -, N 

f q (E) = f q (E + U/2) « — / dE'f q (E') « - £ UE)- (26) 

^ Je + AE/2 IS/ ^[ 



The width of the box average, AE, is related to the width the Lorentzian energy average by AE s=s ttI/2 and to the 
mean spacing, d, of the N ND states by AE « Nd. This approximation is good as long as T/d ^> 1 [see Eq. (|39| ) 
below for the definition of T] . Within these assumptions we see from Eq. (jj^) that the go q and g q o are related to Woo 
by 

— ■ — 2i \ - —ni \ - — 7T2 

Woo = — 2^90q9q0 ~ 2^9o q g q o ~ —go q g q a- (27) 
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Thus the spreading with is given by 



and the energy shift by 



277 

T x = T Re go-q-g^ (28) 



A^^Im^. (29) 



From Eq. (E3) we can calculate the amplitude autocorrelation function 



1 , ZZ ; JK\ ST~^ 90qgq09 0q ,g q/ o -T— TTT^ 

— ^Aoo e + e/2 > g , g /0 F I rf io\ li I /o 5 A o(e-e/2) , 31 

Is , (e + e/2 - £q + «r,j/2) (e - e 2 - E q < - iT q ,/2) 



where we have made the variable changes e — E — E' and e = (E + E')/2. Consider the middle factor in Eq. (|3] 
which we anticipate is a function of e only 



/ \ = goqgqog q'g q >o , , 

W ^(e + e /2-£?,+»r,/2) ^-6/2-^-^/2)- V ; 
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We interpret the energy average in Eq. ( p2[ ) to be an average over e. Employing a box energy average 

r E +AE/2 



1 <._C -t-iiB/Z 



90q9q09 q'9 q ' 



(e + e/2 - £ 9 + iT g /2) (e - e/2 - £ 9 < - iT^/2) ' 



(33) 



If AU is large compared to the T g but small enough for the E q and T g to be treated as constants then we may extend 
the limits of integration to ±oo and perform the integral using the calculus of residues to obtain 



a e 



2iri ^ 



90q9q09 q'9 q '0 



AE^f e + E q , -E q + i{T q ,+T q )/2 



(34) 



Assuming that the phases of the go q and the g q o are randomly distributed as a function of q the double sum in Eq. fl34 
collapses to a single sum giving 



a{e) 



2m ^ |ffo g 5go| 2 



AE 



E 



e + «r 



(35) 



Then employing the definition of the average given by Eq. (pq) 



2tR 

a( £ ) « — 



lffOg5 9 o| 



Assuming that the average of a ratio is equal to the ratio of the averages we get 



a{e) 



2ir e + iT' 



(36) 



(37) 



with and T introduced according to Eq. ([38|) and Eq. ( p9| ) below. The introduction of is based on the assumption 
that 



\go q g q o\ 2 « 2 |30gffgo| 2 ~ 2 



2n 



(38) 



The factor of 2 which appears in the first manipulation in Eq. (^) accounts for the self-correlation (present since the 
entrance and exit channel are both |0)) and is equal to the elastic enhancement factor for compound elastic scattering 
in the overlapping resonance region. In the second manipulation in Eq. (^) we have again ignored the energy shift 
[see Eqsj||) and (§§)]. 
Eq. (p7|) also introduces the correlation width 



r 



r 



= -2lm(q\H QQ + W QQ \q) 
where we used Eq. (|l^) for the electromagnetic width of the |Q) and introduced 



(39) 
(40) 
(41) 
(42) 



rT = -2lm(q\W QQ \q) 
= -21m 



9q090q 



E - E a +i{T s +Ti)/2 

v l d (r s + r^) 

2tt (-B-s ) 2 + (r s + ri) 2 /4' 



(43) 



which is the width for their decay back to |0). The approximations represented by Eqs. ( p9| ) and (|42j) will be discussed 
in Section [v|. Using Eq. (^) and approximation ([42] ) we finally obtain 



AUE)AUE'f « 2(2^r Ar /d)- 1 (rVr s ) 2 a qo {e) 



iT 



N 



E -E' + iT N 



(44) 
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When E' = E this reduces to 



I | = 2 (2 7 rr JV /d)- 1 



r 2 r i 2 



[{E-E y + (T S + Tiy/AY 



(45) 



which is the average of the fluctuation contribution to the transition intensity. 

The integrals in Eqs. ( |l6| ) and ( |l7|) for the average and fluctuation contributions to the total decay intensity may 
now be carried out using the calculus of residues. Substituting Eq. (|2l]) into Eq. (|l6|) we obtain 

1 



J av = (46) 

Eq. fllq) is identical with the equivalent result in GW (see also Ref. [fL4|). Substituting Eq. d45| ) into Eq. (|l^) we 
obtain 



-i (rVr s ) 2 
(i + n/r s ) 3 
= 2(^r Ar /d)- 1 / i a n v (i-0 2 . 

for the average fluctuation contribution to the average decay intensity. 



(47) 
(48) 





FIG. 1: Average of the fluctuation contribution to the intraband intensity vs. log 10 (bj) where bj = /Ts- The solid lines 



were calculated using Eq. ( |47[ ) and the dotted lines by GW's fit formula, Eq. (|49|). The variable IV/d took the value 0.1 for 
graph (a) and 1 for graph (b). 



Eq. @ for J? is plotted in Fig. and for comparison we have also plotted a fit formula which was obtained by 
GW, 



It 



1 - 0.9139 (T N /d) 



0.2172 



exp 



0.4343 In 



(T N /d) 



0.45 



-0.1303 



2 ^ 



-0.1477 



(49) 
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Qualitative agreement is seen between the two formulae. Note that Eq. ( ^9[ ) yields a negative intensity for IV /d > 1.51 
which excludes its use in this regime. Our result which is only strictly valid when Tpj/d 3> 1 is simply inversely 

proportional to T^/d. The exact result of GW for J|* [Eq. (24) in GW] can be used for any Y^/d also decreases 
monotonically with increasing Y^/d. 

The dependence of 1^ (and that of on r^/rs results from the resonant doorway energy dependence of the 
decay amplitude ^4oo (E) [Eq. (pi])] . This energy dependence also manifests itself in the average of the fluctuation 
contribution to the transition intensity | Aq (i?) [ 2 [Eq. (fig)]. GW include precisely the same energy dependence in 
their calculation by use of an energy dep ende nt transmission coefficient to describe decay to the SD band [see Eq. J7S| ) 
below and the discussion in Subsection |V Cf| . This is the reason for our qualitative agreement with GW concerning 



IV. VARIANCE OF THE DECAY INTENSITY 



The error incurred in making the energy average is given by 

A J — T —7~ — — T R 4- 2Re 



dEA 00 (E)A% (E) 



(50) 



The average of the error vanishes: A/i„ = 0. A measure of the dispersion of the calculated Zi n is given by the variance 

r ^ (51) 



(A4 



/in - 40 



In order to evaluate (A/j n ) the averages indicated in Eq. (pl|) must be performed before the integration which appears 



in the definition of /;„, Eq. (|1|). We obtain 

-2 



(A/ in ) z « (27rr s ) 
In deriving Eq. 



dE 



dE' 



A* (E)AUE>)* + 2Re A m (ET ^ A 00 (E>) 



(52) 



we have used A% o {E)=A% o (E)A% o (E')=0, A$ (E)A$ (E') ^0 and 



AUE)AUETAUE')A$ (E') = Ko^Ol + ^Ko^O 

ages of terms containing odd powers of ^4p n vanish. 

Substituting Eqs. (Q) and (|2l]) in Eq. (|5^), making the changes of integration variable 

E-E a = {Y s +Y±)x/2 

e'-e = (r s + rV/2 



and we have assumed that aver- 



and using Eq. (46) and (E^) for I™ and J? we are able to write (AJi n ) in the form 



where the variable £ is defined by 



£ = 



(A40 2 = i& /l (0 + 2«/2 (0 : 



^-(l + rVrs)^/^- 1 

1 Af 



= ^(i + rs/r'r 1 = ^(1-720. 

1 JV Ijv 
The functions /i and /2 which have been introduced in Eq (p5[) are given by 



/i(0 = 



da- 



eta' 



oo [(x - x'f + A/e\ [x 12 + l] 2 



(53) 
(54) 



(55) 

(56) 
(57) 
(58) 

(59) 



7T 2 £ 

-4 



Im 



dx 



Im 



-co [x 2 + 1]' 

dx 
\x 2 + ll 2 



[jj - .t' + 2i/£] [x' + i] [x' - i] 
1 i 



+ i(2/£+l) { x + i(2/Z+l)Y 



(60) 
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and 




log 10 © 

FIG. 2: The functions /i (£) [Eq. (||)] (solid line) and / 2 (£) [Eq. @] (dotted line) plotted vs. log 10 (£). 



and 

The integrations in the calculation of /i and f% above were again carried out using the calculus of residues and were 
checked by numerical integration. The functions f\ (£) and fi (£) are plotted in Fig. ||. 

We have thus shown that a complete description of the decay-out of a superdeformed band within the energy 
average approach requires three requires three dimensionless variables, T^/rs, Tjy/d and r7v/(rs + T^). We find 



(1) The average contribution of the background to the intraband decay intensity, If£, [Eq. (|46|)] depends only on 

rVr s . _ 

(2) The average of the fluctuation contribution, lf n , [Eq. (Ejj)] depends on two variables: T^Ts and Tn jd. 

(3) The variance, (A7 in ) 2 , [Eq. @] depends on three variables:^ /Ts> Ly/d and Fjv/(r s +T l ). 

Fig. H shows a plot of the average intraband decay intensity I[ n [Eq. (|l5|)] calculated using Eqs. ( (46[) and ( |47|) for 



1^ and 7?. For comparison we also show the L ln which results when the GW fit formula [Eq. (49)] is used for 7? 



instead of Eq. (Ed). The two curves are barely distinguishable for T^/d — 1. The GW fit formula incorrectly gives 
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-3-2-10 1 2-3-2-10 1 2-3-2-10 1 2 3 

log 10 (bj) log 10 (bj) log 10 (bj) 



FIG. 3: Average intraband intensity 7i n vs. 

l°Sio( fo j) where bj = T l /F s . The filled circles were calculated using Eq. ( |lq ) 

together with Eqs. ( |4rj| ) and (|47[). The error bars for the filled circles show \J (A/i n ) 2 calculated using Eq. ( p^) with with £ 
in the form given by Eq. Q5q). The dotted lines were calculate in the same manner as the filled circles except that GW's fit 
formula, Eq. (fi9|), was used in the place of Eq. ([j?]). The variable Fjv/d took the following values: 0.1 in graphs (a), (b) and 
(c); 1 in graphs (d), (e) and (f); 10 in graphs (g), (h) and (i). The variable T^/Tn took the following values: 10 -3 in graphs 
(a), (d) and (g); 1 in graphs (b), (e) and (h); 10 3 in graphs (c), (f) and (i). 



intensities which are greater than unity when Tpj/d = 0.1 as does our Eq. (f47|). The exact formula of GW [Eq. (24) 
in GW] does not suffer from this problem. Our results are only strictly valid when T^/d 3> 1. 

The error bars in Fig. || were calculated using Eq. (|35|) with £ in the form given by Eq. (j58[) . We calculate 



(A/jn) 2 for three values of the ratio /Tn (see figure caption). Since (A7 n ) 2 is proportional to 7? it also increases 
monotonically to a maximum before decrease monotonically to zero as a function of T^/r^. For the same reason the 
error increases monotonically with decreasing Tpj/d. The error estimate presented in GW exhibits the same trends 
with r^/rg and Tpf/d. 

Since the variance depends only on (Ts + T^)/Tn in addition to r^/rs and Tpj/d, upon fixing the latter two 
variables the variance may be considered a function of any one of T l /T N , T S /T N , T l /d or T s /d [see Eqs. @ and 
(|58|)]. In Fig. H we chose to use Eq. (|5|) for £ and fixed the value of T^/r^v- When Ts/^n is fixed instead, a slightly 
different dependence on T^/r^ is obtained [compare Eqs. ( |57| ) and (p8|)]. Fig. ^ shows a plot of the standard deviation, 



(A/i n ) , [Eq. (p5|)1 as a function of T^/Tn for fixed r^/rs and TN/d. Ultimately, the variance like the intensity is 
a function of the spin of the decaying nucleus and could provide an additional probe to the spin dependence of the 
barrier separating the SD and ND wells which is contained in the spreading width ||, ^, [?], ||, ^] . 

Our result for the variance of the decay intensity, (A/; n ) 2 [Eq. (55)] has a structure reminiscent of Ericson's 
expression for the variance of the cross section [0 . This connection will be fully explored in the Section For now 
we note only that what distinguishes Eq. ( |55| ) from Ericson's expression for the variance of the cross section are the 
functions fi (£) and fa (£) which result from the energy integrations in Eq. (E3). 
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0.24 r 

0.22 




-5-4-3-2-10 1 2 3 4 

log 10 (Cj) 



FIG. 4: The standard deviation of the decay intensity y (Aim) 2 vs. log 10 (cj) where cj = I^/Tat plotted using Eq. ( |5E| ) with 
£ in the form given by Eq. (^) for fixed bj = r^/Ts and Tjv/d. 



V. VARIANCE OF THE DECAY INTENSITY VERSUS AUTOCORRELATION FUNCTIONS OF 

STATISTICAL NUCLEAR REACTION THEORY 



A. Results familiar from statistical nuclear reaction theory 

All moments of the ^-matrix, S a b(E), the quantities that describe the way in which S a b(E) fluctuates about it's 
average, S a b{E), can be expressed in terms of S a b(E) itself |19] ], Normally, specific moments such as the amplitude 
and cross section autocorrelation functions are expressed in terms of transmission coefficients, defined to be 

T a = T aa , < T a < 1, (65) 

and their generalisation 



T a b = 1 - J2 S «> S >>c*=J2 S %c S V- (66) 

c c 

Here S^ b = S a b — S a b is the fluctuating part of the S'-matrix. The transmission coefficient T a is the probability of 
transmission from a compound nucleus resonance to channel a and is obtained from the optical model. 

In what follows we quote results for both the amplitude and cross section autocorrelation functions in the overlapping 
resonance region for the purpose of comparison with our results for the decay intensity. The amplitude autocorrelation 
function c a b(E,E') for amplitude 

A a b = S ab - iSab, (67) 

is defined by 



Cab{E,E') = Aab(E)Aab(E')* - Aab(E) Aab(E')* (68) 



= A%{E)A%{E'Y (69) 

- ^bZWTv- (70) 
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The correlation width T is given by 



Derivation of Eq. (70) requires the assumption that A ab (E) = A ab (E'). This assumption together with statistical 
assumptions equivalent to those employed in our treatment of the factor a(e) [Eq. (j32[)l gives rise to a c ab {E, E') which 
depends solely on the difference of the two energies E — E' . 

Eq. (HH), which is essentially cqq{E, E') does not depend solely on E — E' . The background energy modulation in 
both E and E' which is characteristic of an isolated doorway resonance is explicit. In the present case it cannot be 
assumed that Aqq(E) ^ Aqq(E') for arbitrary E and E' . The double integral for the variance of the decay intensity 
[Eq (52|) and Eq. ( [77] ) below] is sensitive to this fact this as it contains products of the these background amplitudes 



at ar 



Ditrary E and E' . 



The amplitude autocorrelation function, Eq. (|3|), contains two distinct energy dependences, one characterised by 
r 7v which is analogous to Ericson's correlation width as defined by Eq. (|7l]) and another characterised by Tg + the 
width of the doorway. Writing Eq. (|44| ) in terms of x and x' defined by Eqs. (J5^) and ( |54|) 



16? — 1 

c 00 (E, E') = -^/M g a, (72) 

i m ln [x + i\ 2 [x-x' + 2i/£\[x> -if 

we see that it in fact depends only on the ratio Tjf/(Ts+T^) and it is through Eq. ( [flf ) that this variable enters our 
calculation of the variance of the decay intensity. 

The amplitude autocorrelation function is not an observable quantity. The correlation width, Eq. (|7l|), must be 
extracted from correlation analysis of the cross section. The cross section autocorrelation function, C a b(E, E'), for 
cross section a ab = |^4 a b| 2 is defined by 



C ab {E,E') = a ab {E)a ab (E>) - a ab {E) a ab (E>). (73) 
« \c ab (E, E')\ 2 + 2Re A ab (E)* c ab {E, E')A ab {E') (74) 



7 ah + LO ab a ab 



{E-E'Y+Y 2 ' 



(75) 



Here, cr^ = |^4 a fc| is the background cross section. The fluctuation contribution to the cross section in terms of the 
transmission coefficients is given by the Hauser-Feshbach formula 

_fl .. T a T b (7g) 



or some modification of it designed to account of width fluctuations, direct reactions etc. |T5|, pQ |. 

Eq. ( |52"| ) for the variance of the decay intensity can be written in terms of thecross section autocorrelation function 
defined by Eq. @ as 

/oo poo 
dE / dE'C o{E,E'). (77) 
-00 j — 00 

The same comments concerning the energy independence of the background amplitude apply to the derivation of 
Eq. (j7^) as apply to the derivation of Eq. (]70|). Likewise Cqq{E,E') in the case of the present paper [the integrand 
in Eq^52|)] is distinguished from Eq. ( f75|) by it's explicit inclusion of the energy dependence of the background. 



Eqs. (|52) and ( |74| ) assume that that only pairwise correlations are present. Both Eqs. (f70|) and ([75]) are valid when 
Y] T a 3> 1, that is, in the strongly overlapping resonance region. 

B. Expression of the decay intensity and variance in terms of transmission coefficients 

Following GW we introduce two transmission coefficients, Tq(E) and TV, where 

T (E) = 1 - I = {E _ Eo) 2 T f {rs + Tl)/i - (78) 

4 W-^ V ) (79) 



A(E-E )y(r s + ri) + v 
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describes transmission from the \Q) to the SD band and 



Tv = 2nT N /d, 



(80) 



describes their transmission to ND states of lower spin. We have not derived Eq. (|8C|). For the purposes of the present 
paper it can be taken as the definition of Tv. The reader is referred to the discussion in Section VIII. H of Ref. ]2lf| 
which contrasts the relation of the correlation width, T, to transmission coefficients with the the corresponding relation 
for the average width, T q . 

We have written T (E) in the form given by Eq. ( pZOj ) in order to emphasise that it is not simply a function of a 
single dimensionless variable, the ratio r^/Fg. It is energy dependent, the energy dependence being characterised by 
Ts + T- 1 -, the total width of doorway state |0). Only it's maximum T (E ) = 47-^(1 — 1^) can be expressed solely 
in terms of T^/Ts- Thus, a quantity sensitive to the gross energy dependence of Tq(E) should depend on Ts + 
Writing the average decay intensity Eq. ( |l5|) in terms transmission coefficients 



1 - (2ttI\ 



dE 



{t (T)-2[T (£)] 2 /Tv}, 



(81) 



we see that it compares the total width of the doorway |0) with the width for the feeding of |0) [thanks to inclusion 
of the normalisation factor 27rr,s in the definition of I m in Eq. (|l|)]. The variance, Eq. (|52|), may be written 



/OO />OG f 

dE / dE'l — 
-oo J -oo y [2ir 



AT (E) 2 T (E') 2 



(E-E')/dY 



T 2 

N 



41m 



A QQ (E)T (E)T (E>)A 00 (E>y 
2tt(E- E')/d + iT N 



(82) 



As discussed in Subsection V A and made explicit by Eqs. ( |59|) and (|6l|), the integrand of Eq. (|82| ) which clearly 
contains two characteristic energy scales in fact only depends on their ratio, IV /(r,g + ]>), the ratio of the correlation 
width to the doorway width. 

Eq. ( |4l| ) may also be expressed in terms of the transmission coefficients T (E) and T/v. Using Eq. ( |4^ ) for T^ we 
get 



2tt~ 



-T, 



N 



1 - 



T (E)/T N 



(83) 



so that the neglect of T T in V is justified when T (E)/T N < I? 
Eq. (pq), in terms of the transmission coefficients: 



< 1. Let us also write the correlation length £, 



£ - (2n/T N ) 



r, 



(84) 



In the case compound nucleus scattering, extraction of T from a measurement of cross section autocorrelation func- 
tion, using say Eq. (|75|), permits the determination of the density of compound nucleus states, 1/d, by application of 
Eq. ( |7l| ) p2] | . A more recent example of energy-autocorrelation analysis of may be found in Ref. p3]| where fluctuations 
in dissipative binary heavy ion collisions are studied. In the present case of the decay out of a superdeformed band 
extraction of £ from thevariance of the intensity, permits the determination of the ratio (Ts +T^)/Tn, or, given T/v 
(equivalently T/v/d) determination of the ratio (Ts +T^)/d. 



C. Comparison with the results of Gu and Weidenmiiller 

GW also take inspiration from statistical nuclear reaction theory but use the MPI approach . The MPI approach 
is concerned with the analytic calculation of ensemble averages, a procedure which is equivalent to the calculation 
of energy aver ages. Ref. [p4[ use the supersymmetry method of calculating ensemble averages to derive an exact 

expression for S^ b (E)S^ d (E'). Their result is found to be expressible in terms of the difference of the two energies, 
E — E', and transmission coefficients. The transmission coefficients themselves are expressed as functions of (E+E')/2. 
The relationship between the results of [[15], [l6| , [l 7[ ] and those of is discussed in Refs. [^0[ ^5j. Several results of 
p5| [l6| , [T^ l can be obtained from that of Ref. p4T by expanding in powers of the transmission coefficients or inverse 
powers of the sum of the transmission coefficients [gOf . 

Calculation of the average of the fluctuation intensity requires the energy integral of the average of the product 
of two ^-matrix elements at the same energy. GW use the results of Ref. [pi for \Sq (E)\ 2 to calculate the average 



decay intensity. As was already noted in Section III , GW include the energy dependence of the background amplitude 



13 



characteristic of an isolated doorway resonance in their calculation by using the energy dependent transmission 
coefficient T (E), Eq. @, in their Eq. (24) for i£. The fact that we use the same energy dependence as GW for 
the background is responsible for the agreement we obtain with GW concerning i; n 's dependence on T^/Ts- The 
differences between our results and those of GW for the decay intensity stem from the assumptions we make which 
restrict our results to T^/d 3> 1. 

Calculation of the variance of the intensity requires the 4-point function at two energies integrated over both 

energies, that is, it requires Sq (E)Sq (E)* Sq (E')S$ (E')* integrated over E and E' . Calculation of the 4-point 
function at two energies was carried out using the supersymmetry method in Ref. |26[ . Their result, like that of 
Ref. [ p4| for the 2-point function depends explicitly only on E — E' and the transmission coefficients which are again 
expressed as functions of (E + E')/2. Within the assumption that only pairwise correlations are important, as was 
assumed in Eqs. (|5^) and (|74]), the 2-point function is enough to calculate the variance. Ref. [^0| showed numerically 
that the exact expression of Ref. |24j specialised to the amplitude autocorrelation function confirms the correctness 
of Eq. ( f7C| ) in the region of strongly overlapping resonances. However, unlike Eq. (|7^), the amplitude autocorrelation 
function as given by Eq. (Q) depends on the background amplitude at two different energies, that is, it depends 
on Aqo(E) and Aqq(E'). When E = E' it reduces to Eq. ( f45| ) which can be expressed in terms of the transmission 
coefficients Tq(E) and Tjv. Thus the decay intensity can be expressed in terms of these transmission coefficients as 
was done in Eq. (jsil). The applicability of Ref. [E4j to calculation of the decay intensity owes itself to the fact that 
the decay intensity may be expressed in terms of transmission coefficients. 

E q. (jyj ) cannot be written in terms of Tq([E + E']/2) and the same applies to the variance as is apparent from 
Eq. fl82|). Thus it is not clear whether Ref. |24| serves as a means to obtain results corresponding to Eqs. (|4|) and ( |8^ ) 
which are valid for arbitrary / d. It would be an interesting challenge to derive an expression for the variance which 
could be used for any value of Tn jd since for the regions which have been most frequently studied experimentally [Q , 
the A « 150 and A « 190 regions, F N /d <C 1. 

GW do not use the supersymmetry method to calculate the variance. They instead estimate the variance by 
performing a numerical simulation. The analytic structure of the variance was not investigated in GW and their 
results make no reference to the variable T]^/(Ts + T^). Given the close resemblance of the conclusions about the 
analytic structure of the decay intensity which may be inferred from the exact result of GW and our approximate 
result for T^/d ^ 1 it seems probable that the dependence of the variance on rAr/(Ts + I>) which we have found 
for Tn /d 3> 1 persists for arbitrary Y^/d. 



VI. CONCLUSIONS 



In conclusion, we have derived analytic formulae for the energy average and variance of the intraband decay intensity 
of a superdeformed band in terms of variables which usefully describe the decay-out. The formulae given by Eq. ( |47| ) 

for the fluctuation contribution to the average intensity lf n and by Eq. ([35]) for the variance (Aij n ) 2 were derived by 
making assumptions and approximations which are strictly valid only in the strongly overlapping resonance region, 
Tiy/d^ 1. However, these formulae are seen from Figs. B and |^ to work well when Yj^/d=X and provide a qualitative 
description even when rjv/<i=0.1. This means that Eq. ( pL7| ) and Eq. ( |55|) cannot be applied to the mass 150 and 190 
regions where T^/d ~ 0.001 but they may prove themselves of practical use in other mass regions. In any case our 
results clarify the analytic structure of the results obtained by GW. In particular we have revealed that the variance 
of the decay intensity depends on the correlation length Tn/ (Tg + T^) in addition to the two dimensionless variables 
r^/Ts and Tpj/d on which the average of the decay intensity depends. Measurement of the variance of the decay 
intensity could yield the mean level density of the ND states in analogy with autocorrelation analysis of cross sections. 
This work was partially supported by FAPESP. 
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